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We determine the ground state of the two-dimensional, fully polarized electron gas within the 
Hartree-Fock approximation without imposing any particular symmetries on the solutions. At low 
electronic densities, the Wigner crystal solution is stable, but for higher densities (rs less than ^ 2.7) 
we obtain a ground state of different symmetry: the charge density forms a triangular lattice with 
about 11% more sites than electrons. We argue that this conducting state with broken translational 
symmetry remains the ground state of the high density region in the thermodynamic limit giving 
rise to a metal to insulator transition. 
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The two-dimensional homogeneous electron gas is one 
of the fundamental models in condensed matter physics. 
Despite its simplicity - the system consists of electrons 
interacting through a 1/r-potential to which a uniform 
positive background is added for charge neutrality - the 
phase diagram at zero temperature is nontrivial[TJ [H |3]. 
In general, it is given in terms of the dimensionless pa- 
rameter = 1/ Y^Trna^, where n is the electronic den- 
sity and aB the Bohr radius. At low density (large r^), 
the potential energy dominates over the kinetic energy 
and the system forms a perfect triangular lattice, the 
Wigner crystal (WC), whereas in the high density re- 
gion {ts 0) the kinetic energy favors a uniform Fermi 
gas (FG) phase [1. Already Wigner [4 argued that the 
unpolarized FG is unstable even in the limit ^ 0. 
Later, Overhauser showed the instability of the unpolar- 
ized WC with respect to spin-density waves, even within 
the Hartree-Fock approximation (HF)[5]. It has fur- 
ther been conjectured that the Coulomb potential pre- 
vents any first order transition between the WC and a 
FG[6]. Despite these rather general instability theorems, 
we are not aware of any quantitative calculations estab- 
lishing the true ground state of the electron gas within 
HF[7]. A previous HF study of the two and three dimen- 
sional electron gas with imposed symmetries considered 
only the FG and different crystal structures for the WC 
phase [8], and only recently an unrestricted HF study of 
the unpolarized three-dimensional electron gas was per- 
formed which proposes a more complicated structure of a 
ground state with spin-density waves in the high density 
region [9 . Indeed, establishing the precise HF phase dia- 
gram of the electron gas is a fundamental question, since 
many-body correlation effects can only be quantified with 
respect to the best HF solution. 

In this letter, we consider the two-dimensional elec- 
tron gas, and, for simplification, we concentrate on the 
HF ground state of the completely polarized system. At 
low densities, ^ 2.7, our simulations always lead to a 



WC, but for higher densities and large enough number of 
electrons, TV, the solution is neither a FG nor a WC: the 
density modulation corresponds to a partially occupied 
crystal of different symmetry compared to the insulating 
WC phase (the number of sites is larger than N). We 
refer to this solution as a metallic phase. 

To observe the metallic phase, the number of elec- 
trons has to exceed a threshold ranging from ^ 10 
at Vs ~ 2.7 up to ^ 10^ at Vg ~ 1. For sizes be- 
low the threshold, the solution is the FG [10]. However 
we can prove analytically that the metallic phase has a 
lower energy than the FG down to = in the ther- 
modynamic limit. This proof also indicates that the 
threshold size varies with as e-x^iC st / r s) ^ which may 
explain why this metallic phase has not been observed in 
previous simulations. It is possible to understand heuris- 
tically, why the metallic phase occurs at small r^. One 
expects that the electrostatic energy is always lowered 
by a periodic triangular charge density. On the other 
hand, at small the kinetic energy dominates and favors 
states close to the Fermi surface. The simplest solution 
is obtained by choosing a triangular modulation lattice 
with generators Qi of modulus 2/ci?. This modulation is 
mainly carried by the wave vectors close to the Fermi 
surface. These new HF solutions open a new perspec- 
tive for the qualitative understanding of the experimen- 
tal observed metal to insulator transition [12 and should 
be considered in studies beyond the HF approximation. 

The N-body Hamiltonian, H = + V, contains the ki- 
netic energy K and the 1/r-periodic Coulomb potential 
V where a uniform positive charge background is sub- 
tracted. Within the HF approximation, the search of the 
true ground state of the quantum many-body system is 
reduced to the simpler problem of finding the lowest en- 
ergy states in the subset of the Slater determinants. Let 
^ = V^i A • • • A (/:?Ar be the Slater determinant associated 
with the single particle states {ipi} and £^((/?i, (/?2, • • • , ^at) 
the corresponding energy expectation value. The varia- 
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tion of E with respect to a variation of the single particle 
state 6(pi is then given by 

SE = ^{h^(pi\S(pi) ^^{S(pi\h^(pi) . (1) 

i i 

where h^^ the so-called HF Hamiltonian, is a single par- 
ticle operator depending on the full state ^ (not on the 
particular choice of the (fi^s). Extremal states must sat- 
isfy the following equation 

h^cpi = ^Cij(pj, (2) 

3 

where Cij are the Lagrange coefficients associated with 
the normalization constraint {(fi\(fj) = 6ij. Conversely, 
if ^ = (pi A • • • A (fN is not an extremum, we have 

h^(pi = ^ Cijcpj + Oi (3) 

3 

where the ^^'s satisfy {(Pi\Oj) = 0, Vi, j. Within the steep- 
est descent method one chooses ffist a N x N unitary 
transformation A = {aij) such that one obtains 

{^'i\0'j} = 0, {0',\O'^}<x{ip',\^'^} = 5ij, Vi,j (4) 

for the transformed single particle states 6[ = ^jCiijOj^ 
= Y^j (^ij^j- The energy E{(pi + A6>i, cpN + ><0n) 
can be expressed as a sum of rational fractions whose 
numerators and denominators are polynomials of order 
four, at most. Thus, it is possible to find the best A and 
to iterate the process until a stationary state is reached. 

In fact, this method has the same drawbacks as the 
steepest decent method in linear optimization problems; 
in general, it converges slowly. For linear problems, con- 
jugate gradient methods are preferable [131 El- How- 
ever, since the HF states do not form a linear space, 
the genuine conjugate gradient method does not apply 
here. We have therefore adapted a variant of this method 
to the non-linear case. Let rji be the previous variation 
5(pi, and Oi is obtained by Eq. We then compute 
E{(pi + A^i + /iT^i, • • • ^(fn-^ A^AT + /i^Ar) for six values of 
the pairs {A, fi} in order to approximate £' by a polyno- 
mial of order two in A and /i. Minimizing the polynomial 
with respect to A and /i, we obtain the new changes of 
the single particle states, Scpi^ and the corresponding en- 
ergy change. This process is iterated until the relative 
variation of the energy, SE/E^ is sufficiently small. 

Since we expect the electrons to crystallize on a tri- 
angular lattice at low densities, we have chosen periodic 
conditions compatible with this geometry. Thus, the unit 
cell of our periodized system is given in terms of two vec- 
tors {Li, L2} of length L and with an angle of 60 degrees 
between both; the volume of the unit cell is 1] = I/^\/3/2. 
We have restricted our study to system sizes which are 
compatible both with the triangular lattice and with a 
closed shell occupation in /c-space. Any triangular crys- 
tal with unit-cell vectors {61,62} is compatible with the 
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FIG. 1: Extrapolated energies —^FG (milli-Hartree units) 
versus Vs- Points with error bars: present calculations, full 
line (red): data of Ref. [8], full line (green): fit to present 
results for (rs ^ 2.7), vertical dash line (blue): ~ 2.7. 



boundary conditions if it satisfies Li = lei -\- me2 and 
L2 = — mei + (/+m)e2, where (/, m) are two non-negative 
integers. The number of sites of the lattice is given by 
Nc = det(Li, L2)/det(ei, 62) = + + Im. For a 
commensurate lattice N = Nc. 

We compute the wavefunction on a x Ng grid, the 
fast Fourier transform is used to switch between real 
and reciprocal space [15 . We have systematically checked 
the convergence of the solution with respect to the grid 
size. For the EG ground state, convergence is reached 
once all /c- vectors up to 2k p are represented in the grid 
{Ng ^ A^N/n). At larger r^, in the WC phase, the 
wave functions are essentially gaussians[T6]. The width 
6 of the gaussians scales as ex (rsN)~^^'^. For a cor- 
rect resolution of the gaussians we need L/Ng ex (5, so 
that the number of grid points increases at low densities, 
Ng cx {NvsY^'^. Convergence was reached for Ng = 32 
(resp. 64, 128) for N < 43 (resp. N < 200, N < 500) 
up to = 30. Whenever the number of grid points was 
chosen too small, solutions without any particular sym- 
metries have been obtained. 

We have further studied the influence of the initial 
state on the final solution, choosing different types of 
wavefunction for initialization: a WC state, a converged 
state stored at larger or lower r^, a state initialized with 
random numbers, or a "metallic state" as described below 
(this state has Nc > N maxima in the charge density). 

Typically, during the minimization procedure, the en- 
ergies decrease exponentially; however, the rate of con- 
vergence depends on the initial state. The decrease in 
energy during transitions to a different symmetry is in 
general much smaller than the convergence within the 
same symmetry. We have often seen energy plateaus 
with changes of relative energy ^ 10~^ just before the 
occurrence of a transition to a completely different state. 
For system sizes up to A/" = 151, the minimization is con- 
tinued until a relative precision of 10~^^ is reached, and 
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FIG. 2: Main points of r(q) Eq. Q for N=151. The black 
(resp. white) dots corresponds to = 3.33 (resp. 2.7). The 
radii of the dashed circles are kp and 2kF- The ratio of the 
radii of the solid circles is -^/IST/ISI. The other points in 
T(q) are several order of magnitude less. 



for larger N a relative precision of 10~^ is used. 

We have studied systems with up to 500 electrons at 
densities corresponding to = 1 up to = 30. In Fig- 
ure [l] we report the energies of the obtained HF ground 
state Eoo{rs)^ extrapolated to the thermodynamic limit, 
as a function of r^. In the high density region, for > 3, 
we obtain good agreement with the results of Trail et 
a/. [2, which imposed the triangular WC symmetry. For 
Ts < 3, we find ground states of different symmetries 
than the WC which we describe below. These new states 
have been eventually found down to = 1, at smaller 
Tg, the FG solution is stable for our finite system sizes 
with N < 500. 

We now turn to a more detailed analysis of the different 
phases. It is instructive to follow the evolution of the 
Slater determinants in /c-space as a function of r^. In 
order to quantify the differences to the FG solution, we 
project each plane wave of the FG {ko < kp), on our 
HF solution {(pi}^ 

= 4o+K^ (5) 

^ko = Xl('/^i.'Ako)¥'i (6) 

i 

where is the projection of on the subspace 
spanned by the (pi^s. The deviation of the HF solution 
from the FG can be measured by r(q) as 

T(q) = E kk„(q-ko)|' (7) 

ko<kF,ko7^q 

For the FG, r(q) is 0. Remarkably, for the WC or the 
metallic ground states, r(q) is non zero only for few q- 
vectors (see Fig. [2|. For the WC, r(q) is nonzero on 
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FIG. 3: Charge density p/{p) - 1, with p{x) = J^f^i 
Top: the number of maxima is AT^ = 92 + 52 + 9 X 5 = 151. 
Bottom: the number of maxima is AT^ = 12^ + 1^ + 12 X 1 = 
157. Gray levels corresponds to the same density in both 
figures. Colored lines correspond to Li = lei + me2, where 
the numbers stands for /, m (see text). 



the points corresponding to the triangular reciprocal lat- 
tice: aqi + 6q2, with a, b integers and qi = /Ki — mK2, 
q2 = mKi + (/ + m)K2 (K^, • = 27r6af3) II]. As 
rs decreases from 26.7 down to 3, only the first layer of 
g-vectors corresponding to the triangular lattice remain 
important. In this region, only the hexagon of the 6 low- 
est (7- vectors of the WC remains (see Fig. |2|, and the 
shape of the charge-charge density in /c-space is hexago- 
nal instead of the circular shape at lower densities. This 
wavefunction has therefore a similar character as the hy- 
brid wave functions studied in Ref. However, in con- 
trast to Ref. [3], instead of a phase transition, we find a 
continuous change from many layers of /c-vectors compat- 
ible with the WC at high to only one layer (Figure [2| 
as Ts decreases. At the same time, r(q) in Eq. Q is dom- 
inated by the contribution of the FG vectors ko close to 
Fermi surface. For |ko| ^ kp^ 0^^(k) has essentially only 
one nonzero component at the vector k = ko + q closest 
to the Fermi surface, with q in the reciprocal lattice of 
the WC. 

At 2.7, the system lowers its energy by delocaliz- 
ing the electrons on a denser lattice with more sites than 
electrons (Figure [2] and [3| . This denser lattice is char- 
acterized by integer numbers (I' ^m') different from the 
WC lattice (l^m). For some system sizes the maxima 
of r(q) correspond to various couples of {I'^m'} leading 
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FIG. 4: Nc — N versus the number N of electrons, where Nc 
is the number of maxima in the charge density. Remark that 
Nc never decreases with r^. The hne is from Eq. ([8|. 



to different number of lattice sites, N' = + m''^ + I'm' 
(FigQ. In this case, our system looks like a periodic 
crystal with an incomplete band in contrast to the WC 
solutions of fully occupied bands, studied in Ref. 
Therefore, our translational symmetry breaking HF so- 
lutions below 2.7, have a genuine metallic character. 
However, as approaches zero, the energy gain of this 
metallic crystal compared to the FG gets more and more 
tiny. At the same time, as one can see in Figj4j Nc is 
either constant or increases when Vg decreases (apart for 
one exception N = 301 at = 2.67). For < 1 and 
N < 500, the FG is always the ground state. 

Now we show that the translational symmetry break- 
ing, metallic states should exist at small in the ther- 
modynamic limit. Let us replace a plane wave state k of 
the free FG (||k|| < kp) by a superposition of two plane 
waves with wavevectors k and k + q (||k + q|| > kp)- 
Choosing q on the six-fold star of a triangular lattice 
we certainly obtain a gain in potential energy. The in- 
crease of kinetic energy is minimized if ||k|| ~ kp and 
||k + q|| ~ kp. Then, the number of solutions for k is 
optimal if ||q|| ~ 2A:i?. This solution corresponds to a tri- 
angular lattice of length Lc = 27r/(>/3/cF) in real space 
leading to a unit cell of volume Qc = VSLl/2. Since our 
system is contained in the volume Q = V^L^/2, we will 
obtain Nc = flc/^ lattice sites, or 
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where we have used n = N/Q = /c|./(47r). From Fig. |4] 
we see, that, as decreases, the number of sites of the 
HF ground state indeed approaches the limit of Eq. ( 
A simple model accounting for such a perturbation [ 
(actually far in norm from the FG) gives an energy dif- 
ference from the FG varying as — exp(— a/rs)/rs at the 
limit Ts ^ 0. Our data are reasonably well represented 
by such a law (see Fig. [T]). 



In conclusion, we have computed the ground state en- 
ergy of the two dimensional fully polarized periodic elec- 
tron gas in the HF approximation. We recover the WC 
at large r^. For ^ 2.7, the ground state changes sym- 
metry. Still breaking translational symmetry, the crys- 
tal formed contains more lattice sites than particles and 
should therefore have metallic properties. For our finite 
systems, this solution disappears at small in favor of 
the FG. The same metallic phase has also been numeri- 
cally observed as the HF ground state both in the unpo- 
larized system and in a model where the electrons inter- 
act with a screened Coulomb potential pT]. 

Post HF approaches, like the Quantum Monte Carlo 
method [U [21 |3], have to be used in order to study the 
influence of correlation effects and to establish precisely 
the phase diagram of the electron gas. 
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